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On hard exclusive reactions in the time-like region 

P. KrolP 

a Fachbereich Physik, Universitat Wuppertal 
Gaufistrasse 20, D-42097 Wuppertal, Germany 

The proton form factor, two-photon annihilations into pp as well as exclusive charmonium decays are critically 
examined. It will be argued that the standard perturbative QCD analysis of these reactions fails, i. e. the need for 
additional contributions can convincingly be demonstrated. Possible dynamical mechanisms such as colour-octet 
admixtures to the charmonium states or diquarks inside baryons, will be discussed and compared to the data. 



1. INTRODUCTION 

At large Mandelstam s and large momentum 
transfer the hard scattering approach (HSA) [jj] 
provides a scheme to calculate exclusive pro- 
cesses. Observables are described as convolu- 
tions of hadronic wave functions which embody 
soft non-perturbative physics, and hard scatter- 
ing amplitudes Th to be calculated from pertur- 
bative QCD. In most cases only the contribu- 
tion from the lowest-order pQCD approach in the 
collinear approximation using valence Fock states 
only (termed the standard HSA) has been worked 
out. Applications of the standard HSA to space- 
like exclusive reactions, as for instance the mag- 
netic form factor of the nucleon, the pion form 
factor or Compton scattering off protons revealed 
that the results are only in fair agreement with ex- 
periment if strongly end-point region (where one 
of the quark momentum fractions, x, tends to 
zero) concentrated hadronic wave functions are 
used. As has been pointed out by several au- 
thors (e.g. Q), the results obtained from such 
wave functions are dominated by contributions 
ifrom the end-point regions where perturbative 
QCD cannot readily be applied. Hence, despite 
the agreement with experiment, the predictions 
of the standard HSA are theoretically inconsis- 
tent for such wave functions. It should also be 
stressed that the large momentum transfer be- 
haviour of the helicity-flip controlled Pauli form 
factor of the proton remains unexplained within 
the standard HSA. 

Applications of the HSA to time-like exclu- 



sive processes fail in most cases (e.g. Gm, F n , 
77 — > Pp)- The predictions for the integrated 
77 — > 7T7T cross-section (| cos6>| < 0.6) are in fair 
agreement with the data whereas the predictions 
for the angular distribution fails. Exclusive char- 
monium decays constitute another class of time- 
like reactions. If the end-point concentrated wave 
functions are employed again, the standard HSA 
provides results in fair agreement with the data in 
many cases. Noteworthy are the failures for the 
decays of the r\ c and the XcO into pp. The stan- 
dard HSA predicts zero decay widths for these 
reactions while experimentally the decay widths 
and the branching ratios are of similar magnitude 
as those of the other charmonium decays into pp. 
The reason for this failure is obvious: the per- 
turbative mechanism produces only pp pairs with 
opposite helicities while the quantum numbers of 
the rj c and the \c0 require pairs with the same 
helicities. It should also be noted that in most 
calculations of exclusive charmonium decays |J 
a s values of the order of 0.2 — 0.3 are employed. 
Such values do not match with a s evaluated at 
the charm quark mass, the characteristic scale 
for these decays (a s (m c = 1.5 GeV) = 0.37 in 
one- loop approximation with Aqc d = 200 MeV) . 
Since high powers of a s are involved in charmo- 
nium decays a large factor of uncertainty is hid- 
den in the predictions. 

In this talk I am going to discuss higher Fock 
state corrections to the standard HSA. Constrain- 
ing the pion wave function from the recent pre- 
cise data on the irj transition form factor p| , one 
observes an order-of-magnitude discrepancy be- 
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tween data and HSA predictions for charmonium 
decays into two pions. In Q contributions from 
the ccg Fock state are suggested as the solution 
of this puzzle. In order to cure the failure of the 
standard HSA for reactions involving protons a 
variant of the HSA has been proposed [|| — Q in 
which the proton is viewed as being composed of 
quarks and diquarks. The latter objects consti- 
tute a particle model for higher Fock state con- 
tributions. 

2. THE 7T 7 FORM FACTOR 

The apparent success of the end-point concen- 
trated wave functions, in spite of the theoreti- 
cal inconsistencies, prevented progress in under- 
standing hard exclusive reactions for some time. 
Recently the situation has changed with the ad- 
vent of the CLEO data on the irj transition form 
factor F Tr7 [Q. The leading twist result for that 
form factor, including a s -corrections, reads ^] 



V2 
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Q 2 



(1) 



[1 + ^W, 



The function K has been calculated by Braaten 
||. /,r is the usual pion decay constant (130.7 
MeV) and (x~ 1 ) is the 1/x moment of the pion 
distribution amplitude (DA), <j>, which represents 
the light-cone wave function of the pion inte- 
grated over transverse quark momenta, k^, up 
to a factorization scale, hf , of order Q. The DA 
can be expanded upon Gegenbauer polynomials, 

3/2 

C n , the eigenfunctions of the evolution kernel 
for mesons 

<ftir(x, = (f>As(x) [1+ (2) 

E b m (%%4Y (2,-i)] 

where the asymptotic DA is 6x(l — x). The 
process-independent expansion coefficients B n 
embody the soft physics; they are not calcula- 
ble at present. The 7„ are the anomalous di- 
mensions, no is a typical hadronic scale, actually 
Ha = 0.5 GeV. Any DA evolves into the asymp- 



totic DA for InQ 2 — > oo. Hence, the limiting 
behaviour of the transition form factor is 



^*7T"y * V^2y*7T / Q 



(3) 



which is a parameter-free QCD prediction. As 
comparison with the CLEO data j| reveals, the 
limiting behaviour is approached from below. At 
8 GeV 2 the data only deviate by about 10 — 15% 
from (|§|). The leading twist result (Q) without 
(with) a s -corrections nicely fits the CLEO data 
for 5 2 (Mo) = -0.39 (-0.17) ± 0.05 and B n = 
(n > 4) lioH- I.e. the required DA is nar- 
rower than the asymptotic one in the momen- 
tum transfer region of a few GeV 2 . The fre- 
quently used Chernyak-Zhitnitsky DA [Q, de- 
fined by B 2 (Mo) = 2/3, B n = (n > 4), is in 
clear conflict with the data and should, therefore, 
be discarded. 

Recently a modified HSA has been proposed by 
Botts, Li and Sterman |l3| in which transverse de- 
grees of freedom as well as Sudakov suppressions 
are taken into account. This approach has the ad- 
vantage of strongly suppressed end-point regions. 
Hence, the perturbative contributions can be cal- 
culated self-consistently. Using a Gaussian for the 
kj_-dependence of the pion wave function 



* w (a;,kj_;/i F ) 



4>k(x,hf) 



xiVexp ( -a 2 - 
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' x(l — x) 



where V = 167r 2 a 2 /(x(l — x)) and, for a DA with 
B„ = for n > 4, a, = 1/(^^8(1 + B 2 ) which 
automatically satisfies the tt° — > 77 constraint 
fbl , ond finds perfect agreement with the CLEO 
data for B 2 (Mo) = -0.006 ±0.014 0,0, i.e. the 
asymptotic wave function works very well if the 
modified HSA is used. 

3. PION DECAYS OF CHARMONIUM 

In view of the results for F ni a fresh analysis 
of the decays Xcj - > tttt is in order. Using the 
information on the n wave function obtained from 
the analysis of F^y , one finds the following values 
for the partial widths 



r(xco(2) 



= 0.872 (0.011) keV 



(5) 
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within the standard HSA || . As usual the renor- 
malization and the factorization scales are iden- 
tified in that calculation and put equal to the 
c-quark mass. The parameter describing the 
Xcj state is the derivative R' P {0) of the non- 
relativistic cc wave function at the origin (in coor- 
dinate space) appropriate for the dominant Fock 
state of the Xcj, a cc pair in a colour-singlet 
state with quantum numbers 2S+1 Lj = 3 Pj. 
m c = 1.5 GeV and, of course, B 2 (no) — —0.39 
are chosen as well as R' P {0) = 0.22 GeV 5/2 which 
is consistent with a global fit of charmonium pa- 
rameters 0] as well as with results for charmo- 
nium radii from potential models |l6| ] . 

In H also the modified HSA is used to calculate 
the Xcj — * kit decay widths. Taking B 2 — and 
the other parameters as quoted above, one finds 

r(Xc (2) -> 7T+7T-) - 8.22 (0.41) keV. (6) 

For comparison the experimental values as quoted 
in 0] and reported in a recent paper of the BES 
collaboration jl8) are 

r(Xco -> tt+O = 105 ±30 keV(PDG), 
62.3 ± 17.3 keV (BES), 

r(x C 2 -> 7T+7T-) = 3.8 ± 2.0 keV (PDG), 

3.04 ± 0.73 keV (BES). (7) 

One notes that both the theoretical results, (||) 
and (^), fail by at least an order of magnitude. To 
assess the uncertainties of the theoretical results 
one may vary the parameters, m c , B 2 and Aqcd- 
However, even if the parameters are pushed to 
their extreme values the predicted rates are well 
below data. Thus, one has to conclude that cal- 
culations based on the assumption that the Xcj 
is a pure cc state, are not sufficient to explain the 
observed rates. The necessary corrections would 
have to be larger than the leading terms. A new 
mechanism is therefore called for. 

Recently, the importance of higher Fock states 
in understanding the production and the inclusive 
decays of charmonium has been pointed out [p^[ . 
It is therefore tempting to assume the inclusion of 
contributions from the |ccs( 3 Si)g) Fock state to 
exclusive \cj decays as the solution to the failure 
of the HSA. The usual higher Fock state suppres- 
sion by powers of 1/Q 2 where Q = m c in the 




Figure 1. Representatives of the various groups 
of colour-octet decay graphs. 



present case, does not appear as a simple dimen- 
sional argument reveals; both the contributions 
to the decay amplitude, the colour-singlet and 
the octet one, behave as \jvr? c . In S the colour- 
octet contributions to the exclusive Xcj decays 
are estimated by calculating the hard scattering 
amplitude from the set of Feynman graphs shown 
in Fig. [l] and convoluting it with the asymptotic 
pion wave function. The colour-octet and singlet 
contributions are to be added coherently. The 
Xcj ~ * 7r7r decay widths are given in terms of 
a single non-perturbative parameter k which ap- 
proximately accounts for the soft physics in the 
colour-octet contribution. A fit ot the data 
yields K — 0.16 GeV 2 and the widths 

r(x c0 (2) -> vr+TT-) = 49.85 (3.54) keV. (8) 

Comparison with (^) reveals that the inclusion 
of the colour-octet mechanism brings predictions 
and data in generally good agreement. The value 
found for the parameter k has a reasonable inter- 
pretation in terms of charmonium properties and 
the mean transverse momentum of the quarks in- 
side the pions. Thus it seems that the colour- 
octet mechanism leads to a satisfactorily expla- 
nation of the decay rates of the Xcj m to two pi- 
ons. Of course, that mechanism has to pass more 
tests in exclusive reactions before this issue can 
be considered as being settled. 
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4. REACTIONS INVOLVING PROTONS 

The standard HSA runs into many difficulties 
with these exclusive reactions as mentioned in 
the introduction. In a series of papers 0— 1| 
a variant of the HSA has been proposed in 
which baryons are assumed to be composed of 
quarks and diquarks. A diquark, being a clus- 
ter of two valence quarks and a certain amount 
of glue and sea quarks pairs, is regarded as a 
quasi-elementary constituent. In so far, a quark- 
diquark state represents an (unspecified) super- 
position of higher Fock states 

\B,X> = *f D (x,\L±)\qD > (9) 

= *£g|0M |WW>+ 

In the diquark model spin (S) and spin 1 (V) 
colour-antitriplet diquarks are considered. As- 
suming zero relative orbital angular momentum 
between quark and diquark and taking advantage 
of the collincar approximation, the valence Fock 
state of a proton with hclicity A and momentum 
p can be written in a covariant fashion (omitting 
colour indices) 



\P;p,X) = f s $ s (x)B s u(p,X) 



(10) 



+ ^=<f v (x)B v ( 1 a +p a /m) l5 u(p,X) 

where u is the proton's spinor. The two terms in 
jT(j| ) represent configurations consisting of a quark 
and cither a scalar or a vector diquark, respec- 
tively. The couplings of the diquarks with the 
quarks in a spin-isospin 1/2 baryon lead to the 
flavour functions 



B s = uS[ u ,d] 

By = [uV {u<d} -V2dV {uM ]/V3. 



(11) 



In the diquark model the following DAs have been 
proven to work satisfactorily well in many appli- 
cations §]-§]: 



$s(x) = Nsx(l - x) 3 exp 
$ v (x) = N v x{l - x) 3 (l + 5.8a; - 12. 5x 2 ) 



X 



1-x 

2 



x exp 



x 



-) 



(12) 



These DAs are suitable adaptions of a meson DA 
obtained by transforming the harmonic oscillator 
wave function to the light cone. The constants N 
are fixed through the normalization convention 
(N s = 25.97 and N v = 22.92). The exponentials 
in (|lj) guarantee a strong suppression of the end- 
point regions. The masses appearing in the expo- 
nentials are constituent masses since they enter 
through a rest frame wave function. For u and 
d quarks 350 MeV and for the diquarks 580 MeV 
are appropriate mass values. It is to be stressed 
that the quark and diquark masses only appear 



in the DAs (12); in the hard scattering kinemat- 
ics they are neglected. The transverse size pa- 
rameter b is fixed from the assumption of a Gaus- 
sian transverse momentum dependence of the full 
wave function and the requirement of a value of 
600 MeV for the mean transverse momentum (ac- 
tually b = 0.498 GeV^ 1 ). The constituent masses 
and the transverse size parameter are not consid- 
ered as free parameters since the final results only 
depend on them mildly. 

Diquark-gluon and diquark-photon vertices ap- 
pear in the Feyman graphs contributing to the 
hard scattering amplitude of a given process. Fol- 
lowing standard prescriptions, these vertices are 
defined as 



SgS 
VgV 



ig s t a (pi +p 2 )t_i 
-ig s t a {g a p(pi +P2) M 
-gpfj. [(1 + K )P2 - npi] a 
-6W [(1 + ~ K P2]p } 



(13) 



where g s = \JAkol s is the QCD coupling constant. 
K is the anomalous magnetic moment of the vec- 
tor diquark and t a = A a /2 the Gell-Mann colour 
matrix. For the coupling of photons to diquarks 
one has to replace g s t a by — \J AitaeD where a is 
the fine structure constant and en is the electrical 
charge of the diquark in units of the elementary 
charge. The couplings DgD are supplemented by 
appropriate contact terms required by gauge in- 
variance. 

The composite nature of the diquarks is taken 
into account by phenomenological vertex func- 
tions. Advice for the parameterization of the 3- 
point functions (diquark form factors) is obtained 
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from the requirement that asymptotically the di- 
quark model evolves into the standard HSA. In so 
far the standard HSA and the diquark model do 
not oppose each other, they are not alternatives 
but rather complements. Interpolating smoothly 
between the required asymptotic behaviour and 
the conventional value of 1 at Q 2 = 0, the di- 
quark form factors are actually parametrized as 



4 3) (Q 2 ) = 



Ql 



Qs + Q 2 
Ql 



Qi + Q' 



(14) 



in the space-like region. The asymptotic be- 
haviour of the diquark form factors and the con- 
nection to the hard scattering model is discussed 
in more detail in In accordance with the 

required asymptotic behaviour the n-point func- 
tions for n > 4 are parametrized as 



a s F<i\Q 



(15) 



P\Q 2 ) = 



ay- 



Qv 



Ql + Q 2 



n — 3 



4 3) (Q 2 )- 



The constants as.v & re strength parameters. In- 
deed, since the diquarks in intermediate states are 
rather far off-shell one has to consider the possi- 
bility of diquark excitation and break-up. Both 
these possibilities would likely lead to inelastic re- 
actions. Therefore, these possibilities are not con- 
sidered in the diquark model explicitly but taken 
into account by the strength parameters. Since 
in most cases the contributions from the n-point 
functions for n > 4 only provide small corrections 
to the final results that recipe is sufficiently accu- 
rate. 

The relations ( |H||l5| ) represent effective param- 
eterizations valid at large space-like Q 2 . It is not 
possible to continue these parameterizations to 
the time-like region in a unique way since the ex- 
act dynamics of the diquark system is unknown. 
A suitable continuation to the time-like region is 
defined by the following prescription: Q 2 is re- 
placed by — s in ( |j^ , fi"5|) which ensures the correct 
asymptotic behaviour and, in order to avoid the 
appearance of unphysical poles, the diquark form 
factors are kept constant once their absolute val- 
ues have reached a certain value (c = 1.3) M. 



il 










' — ___J^ G U ( S ) : 


- * I 




4m 2 


Q* G u (Q 2 ) 




• • ■ • -» 

' ' 



5 10 15 20 25 30 35 40 

Q 2 [GeV z ] 

Figure 2. The magnetic form factor of the pro- 
ton in the time-like and space-like (at Q 2 = — s) 
regions. The solid line represents the predictions 
of the diquark model M. The time- like data (o, 
□ ) are taken from [^l],g2), the space-like data (•) 
from |Et|. 



The analysis of electromagnetic nucleon form 
factors constitutes the simplest application of the 
diquark model and the most obvious place to fix 
the various parameters of the model. The Dirac 
and Pauli form factors (the necessary helicity flips 
are provided by the vector diquarks) of the nu- 
cleon are evaluated by convoluting the DAs (02) 
with the respective hard scattering amplitude [pj . 
The parameters are determined from a best fit to 
the data in the space-like region. The following 
set of parameters 



f s = 73.85 MeV, 
f v = 127.7 MeV, 



Ql 
Ql 



3.22 GeV, 
1.50 GeV 2 , 
= 1.39 



as 
ay 



0.15, 
0.05, 



provides a good fit of the data || . The param- 
eters Qs and Qv> controlling the size of the di- 
quarks, are in agreement with the higher-twist ef- 
fects observed in the structure functions of deep 
inelastic lepton-hadron scattering if these effects 
are modelled as lepton-diquark elastic scattering 
pcfl . The predictions for the magnetic form fac- 
tor in both the space-like and the time-like re- 
gions, are compared to the data [^l]-^3| in Fig. 
H Two-photon annihilation into pp pairs has also 
been investigated within the diquark model. The 
prediction for the integrated 77 — > pp cross sec- 
tion is compared to the CLEO data in Fig. ||. 
At large energies the agreement between predic- 
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Figure 3. The integrated 77 — > pp cross section 
(| cos# |< 0.6). The solid line represents the di- 
quark model prediction 0. Data are taken from 
CLEO ||. 

tion and experiment is good. The prediction for 
the angular distribution is in agreement with the 
CLEO data too. 

The diquark model also allows to investigate 
the decay process rj c — > pp. A calculation along 
the same lines as for the other two time-like pro- 
cesses, using the same DAs and the same set of 
parameters, leads to a decay width of 3.88 keV 
H which is in fair agreement with the data [l7| . 
Note that in the pure quark HSA a zero width is 
obtained. 

5. SUMMARY 

The study of hard exclusive reactions is an in- 
teresting and challenging subject. The standard 
HSA, i.e. the valence Fock state contribution in 
collinear approximation to lowest order perturba- 
tive QCD, while asymptotically correct (at least 
for form factors), does not lead to a consistent de- 
scription of the data. In many cases the predicted 
perturbative contribution to particular exclusive 
reactions are much smaller then the data. The 
observed spin effects do not find a comforting ex- 
planation. In some reactions agreement between 
prediction and experiment is found although at 
the expense of dominant contributions from the 
soft end-point regions rendering the perturbative 
analysis inconsistent. 

In view of these observations it seems that 
higher Fock state contributions have to be in- 
cluded in the analysis. However, not much is 



known about them as yet. We are lacking sys- 
tematic investigations of such contributions to ex- 
clusive reactions. A few examples of such contri- 
butions have been discussed in this talk, namely 
the colour octet model for exclusive charmonium 
decays and the diquark model. More work is 
needed. 
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